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Abstract 

Most Particle Induced X-ray Emission (PIXE) data analysis codes are not focused 
on handling multilayered samples. We have developed an open-source library called 
"LibCPIXE", for PIXE data analysis. It is written in standard C and implements 
functions for simulating X-ray yields of PIXE spectra taken from arbitrary samples, 
including multilayered targets. The library is designed to be fast, portable, modular 
and scalable, as well as to facilitate its incorporation into any existing program. 
In order to demonstrate the capabilities of the library, a program called CPIXE 
was developed and used to analyze various real samples involving both bulk and 
layered samples. Just as the library, the CPIXE source code is freely available under 
the General Public License. We demonstrate that it runs both under GNU/Linux 
systems as well as under MS Windows. There is in principle no limitation to port 
it to other platforms. 
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1 Introduction 



The PIXE technique constitutes an excellent tool for quantitative characteri- 
zation of samples. This has been traditionally done with the aid of computer 
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simulation and fitting programs which are able to solve the inverse problem 
of determining the composition of the sample from an iterative fitting of a 
PIXE spectrum. In the currently available codes [1], the fitting routines pose 
limitations on the type of samples that can be analyzed in order to guarantee 
a reasonably non-ambiguous solution. This is the case, e.g. in GUPIX, where 
layered samples can not contain the same element in more than one layer 
[2]. Such a situation is unsatisfactory, specially in fields like materials science 
where many interesting samples require a more flexible definition. 

In order to address this problem, we decided to make use of more advanced 
techniques for inverse problem solving [3] for which we needed a way of sim- 
ulating the X-Ray yields from a given arbitrary sample under general exper- 
imental conditions. In this work we present a collection of computer routines 
(the "LibCPIXE" library) that can be used to perform such simulations. The 
models and databases used in LibCPIXE (based on those from the DATTPIXE 
code [4]) are discussed in the next section and the computational details are 
treated in section 3. After that, some examples are presented of the perfor- 
mance of the library by comparing simulations and selected experimental re- 
sults. 



2 Simulation model 



The primary objective of the LibCPIXE library is to calculate the characteris- 
tic X-Ray emissions produced by the irradiation of a sample with an ion beam. 
The number of detected characteristic X-rays, iV^, emitted by an element, i, 
due to an electronic transition, j, can be written for a layered sample as: 



where the / index accounts for each layer in the sample, Q is the detection solid 
angle, 6j is the efficiency of the detector for the energy of the j emission, N is 
the number of beam particles and Xn is the concentration of the % element in 
the layer I. Tji is the transmittance for the j X-ray energy due to absorbers 
located between the surface of the layer / and the active area of the detector 
while Tij (x) accounts for the transmittance between the position, x, inside the 
/ layer and its surface. Finally, Si is the stopping force for the beam particles 
in the layer /, and is the X-ray production cross section. Note that Si, err- 
and x depend on the energy of the particle, E and that the integral covers the 
range between £j in and E ljOUt , the energies of the particles when entering and 
exiting the I layer, respectively. 



(1) 
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In the current LibCPIXE code, the X-ray production cross sections for protons 
are calculated using the ionization cross sections from the ECPSSR theory [5] 
with semiempirical corrections from Paul [6] and Reis [7] for the K and L sub- 
shells, respectively, in combination with Coster-Kronig fluorescence yields [8] 
and the Scofield K and L line transition ratios [9,10]. Ions other than protons 
are not yet supported but implementing this should not be complicated since 
previsions were taken when coding this first version. 

In order to calculate the transmittances, absorption coefficients for X-rays in 
the energy range from 1 to 30 keV have been obtained from [11]. 

As for the stopping forces (aka "stopping powers"), the ZBL85 values [12] are 
used by default and Bragg's rule is applied for compound targets. However, 
if the default values are not precise enough, any other values (including user- 
defined ones) can be easily incorporated. 

LibCPIXE supports arbitrary samples (including thin, bulk, intermediate and 
layered specimens). In the case of non-thin samples, enhancement of the X- 
ray yield may occur due to photoionization by the X-rays emitted in any 
other part of the sample. The treatment of such effect (the so called secondary 
fluorescence correction) is not trivial in the case of multilayered samples. A 
paper giving all the details about how this is treated in LibCPIXE is under 
preparation and will be published elsewhere. 

The detector efficiency, e, is obtained experimentally by analyzing known sam- 
ples and using the code to determine €j for the main emission lines. The full 
efficiency curve is then obtained by interpolation on these experimental values. 

It must be noted that, currently, LibCPIXE is focused on calculating the yields 
of the relevant lines of a PIXE spectrum but does not try to actually simulate 
a real spectrum and, hence, background is not calculated and no convolutions 
to simulate finite detector resolutions are made. In order to compare the sim- 
ulation results with real data, the experimental spectra can be processed to 
eliminate background, identify the peaks an extract their areas. Work is in 
progress to incorporate into LibCPIXE a tool for automating this process, 
but meanwhile, the publicly available code QXAS [13] can be used to that 
purpose. 



3 Computational details 

The LibCPIXE library, which was originally based on the DATTPIXE v5.3 
code [4], is written in standard C. A Fortran 90 compatibility layer has also 
been developed so that existing Fortran programs can be linked against LibCPIXE 
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too. 



Routines are included in LibCPIXE to facilitate the interfacing: sample and 
experimental parameter descriptions can be read from ASCII files or passed 
directly as function parameters. Similarly, results can be directly returned to 
the caller program and/or be written in ASCII files. The samples and the 
absorbers are defined in a layer-by-layer basis using either atomic or mass 
areal densities for the thicknesses. 



To build a fast library has been a priority and this was achieved by using 
several coding optimization techniques. In the first place, all the databases 
are pre-cached and stopping force interpolation tables are built to save com- 
putation time. The thickness of the internal sub-layers required for numerical 
integration of eq. 1 is not taken to be constant but is dynamically adjusted to 
guarantee small enough variations of Ty , Si and cry . Finally, dynamic memory 
allocation and intensive use of C pointers are done to improve efficiency 



The LibCPIXE library is easily modifiable and scalable thanks to its modular 
design and the fact of being open source (it is freely available under the GPL 
license [14] from http://cpixe.sourceforge.net]). A simple program called 
CPIXE has been developed to test and demonstrate the capabilities of the 
library. CPIXE links against LibCPIXE and has been compiled and tested 
under both GNU/Linux and MS Windows (porting to other systems/platforms 
should not pose difficulties). 



4 Examples with the CPIXE demonstration code 



The CPIXE program has been used to simulate three example cases for which 
experimental spectra were acquired. The spectra used for these examples are: 
a) 1090 MeV H + beam on an pure Au bulk target, b) same beam on a stoi- 
chiometric Fe203 bulk sample and c)1220 MeV H + beam on a Mn4ir thin film 
(562 x 10 15 at/cm 2 ) deposited over a Si substrate. The simulated and experi- 
mental yields are compared in table 1. A further demonstration of the practical 
capabilities of the LibCPIXE (beyond the limits of the CPIXE simple example 
code) is presented in [3]. 

From table 1 it can be noted that the experimental results — deconvoluted 
from the raw spectra using QXAS — are reasonably well reproduced by the 
simulations, specially the K ai 2 and L ai 2 lines. One would typically choose 
these lines for a fitting procedure in the case that LibCPIXE is integrated 
into an automated fitting program (as in [3]). 
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Sample Line Simulation Experiment 





Au-L Q12 
Au-L /3l 


7720 
1880 


7874 
2535 




Au-Lp 2 


1546 


1470 


An 


Au-L 7l 
Au-L^ 


323 
325 


499 
394 




Au-L^ 


155 


186 




Au-L^ 


46 


69 




Au-L 73 


51 


50 


Fe 2 3 


Fe-K 

J. C Il ai 2 


1.04 x 10 5 
15079 


1.09 x 10 5 
14891 




Si - K ai2 
Mn-K ai2 


1.76 x 10 6 
41158 


1.77 x 10 6 
41473 






5929 


6027 


Mn 4 Ir/Si 


li-L ai2 
Ir-L/3! 


2511 
634 


2508 
665 




Ir-i/3 2 


423 


460 




Ir-L 7l 


127 


128 




Ir-L* 


119 


122 




Ir-L^ 


17 


18 



Table 1 

Comparison of simulated and experimental X-ray yields. 



5 Conclusions 



LibCPIXE is, to our knowledge, the first truly free [14] and open-source library 
to be made publicly available to the IBA community for PIXE simulations. 
It is our aim that developers of existing IBA simulation/analysis codes may 
consider to make use of this library to extend their codes into supporting 
the PIXE technique too. LibCPIXE is already functional — as demonstrated 
in this paper and in [3] — but more capabilities and extensions are being 
developed at the moment of writing this. 
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